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ESTIMATING  SPECTRAL  INDICES  FROM  TRANSFORMS  OF 
DISCRETE  REPRESENTATIONS  OF  DENSITY  FUNCTIONS 


1.  INTRODUCTION 

One  of  the  important  needs  of  the  defense  community  Is  to  be  able  to 
evaluate  reliably  the  effect  on  optical  sensors  of  the  disturbed  atmosphere 
resulting  from  a  high  altitude  nuclear  event  (HANE).  If  we  can  understand 
the  structure  and  emission  characteristics  of  the  disturbed  atmosphere  It 
will  be  possible  to  design  detectors  which  avoid  particular  wavelength 
regions  and  can  discriminate  between  targets  and  artlfactual  atmospheric 
phenomena. 

A  great  deal  has  been  learned  In  the  past  several  years  from  the 
research  programs  at  the  Naval  Research  Laboratory  (NRL)  and  elsewhere 
about  the  Instability  mechanisms  that  lead  to  Ionospheric  plasma  structures 
and  the  space-time  characteristics  of  the  resulting  strlatlons.  In  the 
near  future  we  hope  to  clarify  the  cause  of  the  so-called  "freezing” 
phenomenon  for  strlatlons,  make  predictions  about  their  Inner  scale  length, 
etc.  However,  to  use  this  Information  to  design  better  detectors,  we  must 
relate  the  observed  spectral  characteristics  of  the  density  fluctuations  of 
the  emitting  medium  to  the  structures  that  numerical  simulations  and  other 
forms  of  data  (besides  optical)  predict,  and  vice-versa.  This  Is  an 
"Inverse"  problem  and  can  be  Ill-posed  and  yield  nonunique  solutions.  The 
motivation  for  our  studies  Is  to  reduce  nonuniqueness,  etc.  through  a 
detailed  consideration  of  the  procedures  used  In  relating  observed  spectral 
quantities  to  strlatlon  properties. 

We  will  address  several  aspects  of  this  problem  In  a  series  of  papers 

that  are  currently  In  preparation.  In  this  paper  we  consider  idealized 

"scan”  functions,  namely  functions  of  one  variable  that  can  arise  from 
Manuscript  approTcd  January  24,  1984. 


three-dimensional  optical  sources  that  are  observed  by  remote  sensors.  We 
will  establish  the  relationship  between  continuity  properties  of  these 
functions  or  their  derivatives  and  Che  spectral  Index  of  their 
transforms.  We  will  show  how  Inadequate  resolution,  whether  In  measured 
data  or  numerical  simulations.  Introduces  errors  In  estimates  of  spectral 
Indices.  Specifically,  we  will  discuss  the  errors  that  arise  from 
"aliasing"  and  the  number  of  data  samples  required  to  obtain  an  adequate 
separation  of  scales.  We  will  illustrate  these  effects  with  numerical 
examples . 

In  a  second  paper  In  preparation  we  relate  spectral  properties  of 
multidimensional  emitting  sources  to  spectral  properties  of  the  scan 
functions.  Sources  of  constant  emission  Intensity  and  sources  with  finite 
gradients  viewed  from  different  directions  will  be  considered. 

In  a  third  paper,  a  simple  model  of  a  realistic  nonaxlsymmetrlc 
Ionospheric  scriation  Is  constructed,  which  incorporates  the  properties  of 
emitting  structures  established  in  the  NRL  research  programs.  Using  this 
model,  we  examine  the  variation  in  spectral  properties  that  would  be 
observed  by  scans  obtained  from  different  viewing  directions.  We  also 
Investigate  the  sensitivity  of  the  spectral  Index  to  variations  In  model 
parameters.  Subsequent  papers  will  Investigate  multiple  strlatlon  effects 
and  other  properties  needed  to  further  clarify  the  relationship  between 
emitting  sources  and  measured  spectral  Indices. 

A  preliminary  Investigation  of  these  topics  has  been  given  by  Wortman 
and  Kllb  [1].  They  have  added  the  additional  feature  of  a  self-similar  and 
probabilistic  distribution  of  scale  sizes  In  the  density  function  and  a 
thorough  cdmparlson  with  available  data.  However,  they  do  not  focus  on 
uncertainties  In  the  spectral  Indices  resulting  from  Inadequate  resolution 


2 


NOTATION  AND  TRANSFORMS 


2.1  Notation 

We  define  the  direct  and  Inverse  Fourier  transforms 


f  -?f  and  f  f. 


Three  forms  of  f  are  considered.  The  continuous  function  on  x  «]; 

the  periodic  function  on  x  e[-L,  L];  and  the  sampled  periodic  function  on 
X  e[-L,  L],  with  2N  samples  and  spacing  h.  These  and. the  relevant  forms  of 
Parseval's  relation  are: 


f(k)  -  /  f(x)  e  ^^*dx, 


X  e  +»], 


f(x)  -  (2v)“^  /  f(k)  e‘^^*“dx, 


+lkx. 


/  Ifl^  dk  -  2tt  /  f^(x)  dx; 


-Ik  X 


f  =  f(k^)  -  (21)"^  /  f(x)  e  dx,  X  e[-L,  L] ; 

-L 


+■  ^  +lk  X 

f(x)  -  I  f(k^]  e  ''  , 


V  e[— ,  -Ho] 


I  |f^|^-(2L)"^/  f^(x)  dx 


n-N 


-inhk 


F  =  F(k J  -  (2N)  I  f(nh)  e  ",  n  e[-N+l,  N] ; 
^  n— N+1 


(2.1a) 

(2.1b) 

(2.1c) 

(2.2a) 

(2.2b) 

(2.2c) 

(2.3a) 
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2«2  Transforms  of  Elementary  Continuous  Functions  of  Compact  Support 
One  of  our  goals  Is  to  characterize  physical  space  functions  by 


spectral  Indices.  That  Is,  the  power  density  Ifl^  (or  IfI^,  etc.)  has  an 
envelope  that  may  be  characterized  by  k  In  various  regions  j  and  we 
seek  to  define  these  regions  and  find  accurate  estimates  of  pj.  For 
convenience  In  Illustrating  separation  of  scales  and  asymptotic  properties, 
we  will  use  a  function  f(z)  that  is  composed  of  piecewlse-constant,  linear, 
quadratic,  etc.  functions  and  fractional  powers  of  these  functions.  We 
believe  these  functions  are  sufficiently  general  to  Include  the  essential 
features  of  real  scans.  The  spectral  Index  p  will  be  determined  by  the 
particular  functions  that  are  chosen  and  by  the  manner  in  which  they 
Intersect. 

First,  consider  the  transform  of  on  x  e[-“,  «] 


^(3  -  /  (3  e"^^dx  -  (ik)  f(k). 


(2.4) 


In  general,  if  3^f  vanishes  sufficiently  rapidly,  then 


(ajf)  -  (ik)'^  f(k).  (2.5) 

In  the  sense  of  generalized  functions,  the  derivative  of  a  Heaviside  step 
function 


H(x) 


rl  X  >  0 
U  X  <  0  ' 


Is  the  delta  function  5(x).  Their  transforms  are  related  by  4^6 (X))  -  1 

A 

*  (ik)  H.  Thus,  for  functions  of  compact  support  that  are  composed  of 


piecewise  polynomials ,  one  differentiates  a  sufficient  niunber  of  times 
until  delta  and  more  singular  generalized  functions  are  obtained.  (For 
example,  for  a  trapezoid  one  differentiates  twice).  To  obtain  the 
transform  of  f,  we  divide  the  transforms  of  these  singular  functions  with 
singularities  located  at  x^j  by  the  appropriate  power  of  (ik)  (the  number 
of  differentiations)  and  combine  with  appropriate  phase  shift  factors,  exp- 
(ikXgj).  Some  typical  results  are  given  in  Table  I. 

Table  1  —  Fourier  transforms  of  elementary  figures 
p  I  Ibransform 


(2) 

2Ak“^ 

-IZ4. 

e  ^  sin 

z_ 

(I.l) 

(2) 

4Ak“^ 

cos  z^  sin  z_ 

(1.2) 

(2) 

4Ak“^ 

sin  Z4.  sin  z_ 

(1.3) 

(2) 

iA[-k' 

.1  ^-ika  ^ 

2A“^  k“^  e  sin  z_] 

(1.4) 

(2) 

-2Ak- 

^  sin  ka  ♦ 

4AA~^  k”^  sin  r+  sin  z 

-  (1.5) 

(4) 

8AA-' 

k-2  e'“^ 

sin*(z_/2) 

(1.6) 

(4) 

16AA- 

1  k-"  cos 

Z4.  sin^(z_/2) 

(1.7) 

(4) 

il6AA 

k-2  .in 

z+  sin*(z_/2) 

(1.8) 

The  symmetrical  trapezoid  of  amplitude  A  Is  obtained  by  combining 


(1.5)  with  (1.2)  where  in  the  latter,  a  ■  0  and  b  •  a.  Thus 


f  ■  A(b  +  a)(sin  sin  z_)/(z^z_), 


(2.6) 


where 


z^  -  i  Vt(b  +  a).  (2.7) 

The  half-width  of  the  trapezoid  may  be  defined  as  ^-^12  *  (b  +  a)/2,  and  the 
first  null  of  (2.6)  is  at  also  the  Interval  between  nulls 

associated  with  sin  z^.  In  a  similar  fashion.  If  f  is 

f(x)  ■  1,  |xl  <  a, 

f(x)  -  [1  -  (1x1  -  aj^M^]^,  a  <  |xl  <  b,  (2.8) 

f(x)  -  0,  Ixl  >  b, 

where  A  ■  (b  -  a),  t  >  -1.  Note,  f  is  singular  at  |x|  ■  b  if  -1  <  t  <  0 
and  f  has  a  singular  slope  at  lx|  ■bifO<T<l.  The  transform  is  given 
in  Appendix  A  and  for  t  ■  1/2  (an  elliptical  arc)  it  is 

f(k)  -  k"^  sin  ka{2  -  v  Hj^(kA)}  +  it  k"^  J^(kA),  (2.9) 

where  is  a  Struve  function,  as  discussed,  for  example,  in  Reference  3, 


and  is  the  Bessel  function  of  the  first  kind.  If  kA  »  1, 


where  Yj^  Is  the  B«>3sel  function  of  the  second  kind.  Thus 


lim  f(k)  -  (2ir)^^^A(kA)"^^^[cos  (kA  -  |l)  -  sin  (kA  -  sin  kal  . 


kA  »  1 


(2.10) 


Further  properties  of  these  functions  are  discussed  below. 


3.  ASYMPTOTIC  SPECTRA  AND  SEPARATION  OF  SCALES 

3.1  Functions  on  the  Infinite  Line 

In  this  section,  quantitative  results  for  special  functions  are 
presented  and  general  rules  are  induced.  Particularly,  that  the  asymptotic 
spectrum  of  continuous  functions  is  determined  by  the  physical  space 
regions  where  slopes  of  f  change  in  the  most  singular  manner. 

The  spectral  representation  of  variables ,  1 f I  ,  in  nonlinear 
dynamical  processes  often  can  be  represented  by  a  power-law  function  k 
in  the  region  <  k  <  k^,  where  pj  is  called  the  spectral  index  in 

region  j.  (This  excludes  the  dissipative  range  where  some  exponential 
variation  with  k  usually  occurs).  These  representations  are  often  obtained 
by  least-squares  fitting  procedures  which  suppress  oscillatory  effects. 
For  example,  the  envelope  of  maxima  of  piecewise  polynomials  discussed 
previously  can  be  fitted  by  If^l  o  k  when  length  scales  are  sufficiently 

"separated".  For  k  •*■  «,  p  is  called  the  asymptotic  spectral  index. 

2  -4 

For  the  trapezoid  (2.6),  the  asymptotic  spectrum  is  (4/A)  k  ,  or  p  ■ 
4  and  for  (2.8),  the  asymptotic  spectral  index  is  obtained  from  Appendix  A 


p  -  2(t  +  1),  T  >  -  1. 


(3.1) 


For  T  «  1/2  (see  Equation  (2.10))  p  *  3  which  is  intermediate  between  2  and 
4,  the  values  for  the  top  hat  and  the  trapezoid,  respectively.  For  -  1  < 

T  <  0,  f  is  singular  at  x  ■  b  and  0  <  p  <  2  (that  is,  shallower  than  a  top 
hat).  For  0  <  t  <  1,  f'(x)  is  singular  at  x  =  b  and  for  t  ^  1,  f'(x)  is 
everywhere  continuous  and  vanishes  at  x  *  b.  From  these  results,  we  induce 
a  rule  for  piecewise  polynomials  and  powers  of  piecewise  polynomials:  the 
as3mptotic  spectral  index  is  determined  by  the  nature  of  the  singularities 
of  3*^f/3x*^,  where  these  piecewise  polynomials  intersect.  For  example,  if 
the  first  derivative  (q  ■  1)  has  a  singular  slope  somewhere,  then 

2  <  p  <  4; 

if  q  *  1  does  not  have  a  singular  slope  but  q  *  2  does  have  a  singularity, 
then 

4  <  p  <  6;  etc. 

Let  us  now  consider  the  separation  of  scales,  namely  different 
spectral  indices  pj  in  different  regions  of  k,  <  k  t  k^.  For  the 

trapezoid  (2.6)  we  have  three  regions: 

f^(k)  -  (b  +  a),  k  «  k^, 

f^(k)  -  2k"\  k^  «  k  «  k^,  (3.2) 

f.(k)  »  (b  -  a)"^  k“^,  k,  «  k. 


where  ■  2Tr/(b  +  a)  and  k2  •  2w/(b  -  a).  Thus  if  (b  +  a)/(b  -  a)  »  1  we 
have  a  good  separation  of  scales.  For  the  trapezoid-plus-top  hat  with  u  « 
1  shown  in  Figure  1, 


f  ■  (1  -  u)(b  +  a)(sln  z^  sin  z_)/(z^z_)  +  2ub(sin  kb)/(kb),  (3.3) 

and  we  have  four  regions 

%  -  (1  -  U)(b  +  a),  k  «  k^, 

fg  -  2k"^,  kj^  «  k  «  k2, 

(3.4) 

fg  -  (b  -  a)“^  k“^,  k^  «  k  «  kj, 

«  2uk"^,  kj  «  k; 

where 


kj^  *  2Tr/(b  +  a),  k2  ■  211 /(b  -  a), 


k^  -  (1  -  u)  li  ^  Tr/(b  -  a). 


(3.5) 


and  where  we  have  assumed  u  «  1.  The  asymptotic  spectral  index  is  2 
because  of  the  small  but  finite  jump.  Note  that  the  last  region  begins  at 
a  point  dependent  on  the  size  of  the  discontinuity,  which  in  practice  could 
be  related  to  a  data  artifact. 


3.2  Periodic  Continuous  and  Discrete  Functions 


Because  of  the  computational  efficiency  of  the  fast  Fourier  transform 
algorithm,  one  usually  imbeds  functions  in  a  periodic  domain,  -L  <  x  <  L. 
If  the  functions  have  compact  support  over  a  range  <  2L,  then  from  (2.2a) 
the  transforms  have  the  same  form  but  the  continuous  k.  is  replaced  by  > 
irv/L,  where  v  takes  on  all  positive  and  negative  integers.  As  a  rule,  if 
one  wishes  many  harmonics  between  nulls,  one  requires  where 

£^^2  "half-width"  of  the  function.  If  we  satisfy  this  criterion  we 

will  obtain  a  reasonable  approximation  to  the  continuous  transform  function 
but  it  may  not  yield  a  good  estimate  for  the  spectral  index,  as  we  will  see 
in  Section  3.3. 

For  discrete  functions,  defined  at  intervals  h  =  L/N,  Eq.  (2.3a)  is 
applicable.  The  discrete  system  has  2N  independent  harmonics  v  »  (-?H-1), 
(-JH-2).. .-1,0,1, ...N.  The  lowest  harmonic  is  (x/Nh)  and  the  highest 
is  (x/h).  For  the  symmetric  trapezoid,  Eq.  (2.3a)  yields 


F(0) 


(B-Hi) 

2N 


sin  9(8+a) 
(S+a)sln(i  9) 


sin  9(B-o) 
(B-a)sin(Y  9) 


•} 


(3.6) 


where 


b  ■  Bh  and  a  ■  oh. 


(3.7a) 


9  ■  kh  ■  xv/N,  V  e[-JH-l,N], 


(3.7b) 


and  9  •  x.  The  essential  difference  between  (3.7)  and  (2.6)  is  the 

flldX 

2  2 

presence  in  the  denominator  of  (sin  9/2)  instead  of  (9/2)  .  This 
difference  is  called  "aliasing"  [4]  and  is  the  result  of  "folding"  the 
discrete  spectnus  of  the  Fourier  series  around  the  highest  mode.  Thus 


aliasing,  an  unavoidable  result  of  dealing  with  discrete  data,  modifies  the 
as]naptotic  spectral  Index. 


To  obtain  a  quantitative  measure  of  the  error  we  define  a  ratio  of 
local  "indices"  and  subtract  one,  or 

d[lnlf^|^]/d(ln  9) 

dflnlF  l^l/d(ln  9) 

L  e 

where  f^  and  are  the  envelope  functions  corresponding  to  f  and  F, 

A 

respectively.  Thus  If  e  (9)  >  0,  F  has  a  smaller  effective  p  than  does 

P  e 

A 

fg.  An  approximate  result  for  the  trapezoid  Is  obtained  by  setting 


1} 


(3.8) 


c  (9)  =  { 


where 


P2  -  2  for  2x/(b+a)  «  k  «  2x/(b-a) 


Pj  ■  4  for  k  »  2x/(b-a). 


Thus 


Sp  -  [(tan  j  9]/(j  9)]  -  1  -  (9^/12)  +  o(9^].  (3.9) 

Is  positive  and  Independent  of  pj  and  Is:  0.024  at  9  *  7r/6;  0.055  at 
9  ■  «/4;  0.103  at  9  ■  it/3;  and  0.273  at  tt/2.  That  is,  aliasing  errors 

decrease  the  measured  spectral  Index.  Thus,  If  we  use  a  nonlocal  fitting 


procedure  to  estimate  p  (as  described  below)  and  we  wish  to  avoid  using 
data  that  contributes  local  errors  >  27Z  (or  >  10.3%),  we  must  discard  half 
(or  two-thirds)  of  the  modes! 

3.3  Estimating  Spectral  Indices  of  the  Trapezoid 

As  discussed  in  Sec.  4,  a  least-squares  (nonlocal)  fitting  procedure 
is  used  to  estimate  spectral  indices.  The  essential  caveats  are:  avoid 
using  data  near  a  transition  between  spectral  ranges;  and  discard  data 
above  (or  .  We  discuss  the  choice  of  appropriate  data 

fitting  regions  for  the  trapezoid,  if  we  wish  to  obtain  estimates  of  P2  *  2 
and  P2  *  4. 

We  wish  to  fit  peaks  of  the  slow  oscillation  associated  with  sin 
k(b-a)/2  in  the  last  region  (No.  3).  The  last  data  mode  will  be 

where,  for  example,  to  avoid  aliasing  errors  y  <  1/2.  If  we  start  at 
3  k2/2,  the  condition  for  peaks  beyond  the  transition  yields  a  range 
condition 

(j  +  U3)(2Tr/(b-a))  -  yw/h 


or 


(b  -  a)  -  h(3  +  2y3)/Y. 


(3.11) 


r:  >■!  ng.i  r  p.'g. j.  jj  a  a  ■■.  *:v <■. 


The  intermediate  fitting  region  (No.  2)  starts  after  kj^,  the  first 
null,  and  proceeds  to  k]^*.  Here  k^*  is  chosen  according  to  the  error  made 
as  one  approaches  k2,  that  is  according  to  the  departure  of  [sin  (k(b  - 
a)/2)/(k(b  “  a)/2)]^  from  unity  as  given  in 

[sin  j  k^*(b  -  a)/  j  kj^*(b-a)]^  »  1  -  +  0(£j^^) , 


or 

-  (3)“^^^  k^*(b  -  a)/2,  (3.12) 

Ne  proceed  as  in  the  highest  range,  and  require  that  we  fit  peaks 
associated  with  sin  k(b+a)/2,  or 

(|+  W2](2ir/(b  +  a))  -  (12)^^^  e^/(b  -  a),  (3.13) 

or  using  (3.11) 

(b  +  a)  -  h(3  +  2p2](3  +  2U3]  ir/Ye^(12)^^^.  (3.14) 

Finally,  we  wish  to  have  sufficient  data  in  the  first  region  before 
the  first  null  at  2N/(S+oi).  Thus 

N  ■  Pj^(S  +  a)/2,  (3.15) 

where  a  minimal  requirement  is  >  4.  If  we  take  y^  ~  4  and  requirements 
in  the  other  regions,  as  follows: 

y2  ■  W3  ■  3.  Y  ■  1/2  and  -  0.1,  (3.16) 
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^  ^  . . 


then  we  obtain 


(b  -  a)/h  -  18,  (b  +  a)/h  -  464.6  and  N  -  929, 

where  the  last  is  obtained  from  (3.15).  In  Section  4.2  we  will  compare 
methods  of  fitting  the  data  from  the  trapezoid  (b-a)/h  *  16  and  (b-l-a)/h  ~ 
464  for  N  «  512,  1024,  and  2048,  which  straddle  the  value  N  >  929. 

4.  FITTING  DISCRETE  DATA 

4.1  Least  Squares  Fits 

In  this  section  we  Illustrate  errors  In  a  weighted  least-squares  fit 
of  discrete  data  lF(k^)l  over  specified  ranges  with 

If  r  -  F^  k"P  -  v”P,  V  e[-  Nfl,  N] .  (4.1) 

Ad  vJ 

Thus,  p  and  Eq  are  obtained  from  the  pair  of  linear  equations 

I  w^  log^glFl^  -  logj^Q  Eg  I  w^  +  p  ^  w^z^  -  0  (4.2) 

V  V  V 

^9  9 

y  w  z  log,  nlFr  -  log,-Ert  Iwz  +pywz^»0  (4.3) 

^  v  V  10  ®10  O^vv^vv 

V  V  V 

where  z^  •  log^gV  and  the  weighting  w^  that  yields  good  fits  Is 

V 

because  it  emphasizes  the  lower  modes. 
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•  i 


The  fitting  range  for  spectral  region  j  Is  defined  as 


V  e 


> 


where  v.^  and  v.„  are  the  Initial  and  final  mode  values  Included  In  the  fit 
ji  jr 

^  2 

for  region  j,  and  are  chosen  where  iF(k^)  I  has  a  local  maximum,  and  such 
that  they  are  not  too  close  to  transition  points.  This  procedure  was  found 
to  give  a  good  estimate  of  p  for  single  figures.  No  consistent 

Improvements  were  obtained  when  spectra  for  single  figures  were  smoothed. 
However,  when  many  figures  were  placed  on  a  line  (Including  overlapping 
figures),  we  found  that  the  precise  Vj^  and  Vjp  were  less  critical.  This 
follows  because  the  point-to-point  variation  of  |F|  was  large  (l.e., 
poorly  correlated).  It  Is  possible  that  an  algorithm  that 

fits  logj^glFl  with  a  polynomial  in  would  give  a  better  estimate  of  a 
local  spectral  Index. 

Table  2  contains  summary  Information  on  the  top  hat  (a  »  b),  circle 
(Eq.  2.8  with  a  >  0  and  t  ■  1/2)  and  trapezoid*  This  table  Illustrates  the 
errors  made  In  obtaining  p,  when  parameters  are  varied  Including  the 
fitting  Interval  [Vj;^  ^jF^  'Hnax*  total  number  of  data  maxima  In  the 

Interval.  Measures  of  the  quality  of  fit  are  given  by  6^  -  (p  -  Pj)/Pj  x 
100  and  by 


a 


[I 


(4.5) 


where  z  ■  ^  w^.  Note:  In  statistical  fitting  procedures  a  Is  the 

V  V 

square  root  of  the  variance  of  p. 
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4.2  Discussion  of  Fitted  Results. 

The  results  given  in  Table  2  are  for  the  top  hat  (cases  l-*6),  the 
circle  (cases  7-LO)  and  the  trapezoid  (cases  11-21).  To  assure  a  good 
value  of  p  (to  minimize  the  figures  were  well  resolved  and  occupied  a 
small  extent  of  the  total  interval  (to  yield  well-defined  oscillations). 
These  requirements  dictate  a  large  total  mesh  size.  The  cases  Illustrated 
in  the  table  show  the  effect  of  satisfying  these  criteria  to  a  greater  or 
lesser  degree. 

If  we  compare  the  errors  in  the  least  squares  results  for  region  2  we 
find  that  for  the  top  hat  and  circle  (two-region  functions)  ^  is  much 
closer  to  P2  than  for  the  trapezoid  (a  three-region  function) ,  for  a  given 
mesh  size.  Thus,  even  with  large-mesh  sizes  a  least  squares  fit  to  the 
data  produces  substantial  errors  in  the  prediction  of  a  spectral  index  if 
there  are  several  distinct  regions  of  k-space  with  different  spectral 
indices . 

Cases  I  and  2  show  results  for  a  mesh  (2N  ■  4096)  much  larger  than  is 
being  used  presently  in  numerical  simulations.  A  plot  of  log^Q|F|  vs. 
log^QV  for  these  cases  is  shown  in  Figure  2.  The  ordinate  is  normalized  so 
that  the  area  under  the  profile  is  unity.  The  lower  scale  on  the  abscissa 
is  defined  in  such  a  way  as  to  provide  a  measure  of  the  size  of  the  object 
relative  to  the  mesh.  The  number  given  at  the  origin  is  the  ratio  of  the 
mesh  size  to  the  size  of  the  figure  (here  it  is  8.8  »  2048/232).  The  first 
null  (which  is  not,  in  general,  a  point  of  the  plot)  occurs  at  1  on  this 
scale,  the  second  at  1/2,  the  third  at  1/3,  etc.  In  Table  2  we  see 
errors,  5^,  in  p  of  3%  and  5%,  respectively,  and  a  consistent  variation 
in  0.  Note  that  Che  fitting  procedure  always  yields  a  p  >  p.  This  is  due 


Co  our  algorithm  for  choosing  the  fitting  regions.  We  are  not  certain 


whether  the  Increase  In  error  from  case  1  to  case  2  Is  related  to  the 


decrease  in  sample  size  (14-1038  to  14-526)  or  has  to  do  with  V2p  in  case  2 
being  farther  from  the  aliasing  region,  since  aliasing  tends  to  spuriously 
reduce  p.  In  any  event,  the  errors  are  acceptably  small.  In  cases  3  and  4 
the  same  top  hat  is  placed  on  a  mesh  half  the  size  of  that  in  cases  1  and 
2.  We  obtain  slightly  better  results,  with  errors  of  2%  and  4.5% 
respectively.  The  minor  improvement  is  due  to  the  particular  choices 
of  V2J  and  V2p  determined  by  our  algorithms.  In  fact,  choosing  in  cases 
1  and  2  to  be  15  Instead  of  14  improves  the  fit  enough  that  cases  1  and  2 
then  show  smaller  errors  than  cases  3  and  4,  respectively. 

With  a  smaller  top  hat  in  a  smaller  box  (2N  -  512,  cases  5  and  6, 
Figure  3)  the  errors  are  larger,  4%  for  case  5  and  11.5%  for  case  6.  The 
normalized  standard  deviation,  o,  has  also  increased.  For  the  circular  arc 
the  results  are  consistent  with  the  top  hat  results.  For  example,  cases  7 
and  8  give  6^  of  2%  and  5%,  respectively.  The  errors  for  cases  9  and  10 
compare  closely,  also,  with  the  analogous  top  hat  cases  5  and  6. 

For  the  trapezoid  (see  Figure  4),  the  errors  in  p2  range  from  17.5%  to 
29%.  Generally,  as  the  data  set  increases,  6^  and  "a  decrease  if  we  do  not 
approach  too  close  to  transitional  points  or  aliasing  regions.  The  smaller 
errors  in  region  3  are  the  result  of  aliasing  errors  competing  with  the 
errors  Introduced  by  the  fitting  procedure.  The  variation  In  Sp  for 
different  meshes  among  comparable  regions  (e.g.  26.5%,  23.5%  and  17.5%)  is 
an  '  indication  of  the  magnitude  of  the  variability  obtained  by  such 
procedures  when  data  with  aliasing  errors  are  included.  These  larger 
errors  result  because  the  configuration  of  "minimal"  parameters  (3.16)  does 


aot  yield  a  sufficient  number  of  oscillations  of  data  far  from  the 
transitional  points*  If  we  choose  a  more  conservative  set  of  parameters, 

e.g. 

1  2 

■  8,  U2  *  ^  *  T  ^1  "  0*01 

then  we  obtain 

(b  -  a)/h  -  33,  (b  +  a)/h  -  3290  and  N  -  13,170. 

This  large  value  of  N  cannot  be  used  conveniently  or  routinely  with 
present-day  computers. 

In  addition  to  the  trapezoid  runs  in  Table  2  we  have  made  runs  with 
randomly  placed  multiple  (three)  identical  trapezoids.  In  a  mesh  with  N  • 
2048,  analogous  to  the  single  trapezoid  cases  11-14,  we  found  only  small 
differences  In  the  resulting  values  of  p^.  This  suggests  that  single 
figures  can  provide  a  good  estimate  of  variability  due  to  fitting 
procedures  in  a  situation  where  multiple  figures  occur.  However,  a  proper 
statistical  theory  for  multiple  figures  is  needed  to  generalize  our  limited 
findings. 
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CONCLUSIONS 


We  have  demonstrated  the  essential  elements  which  are  required  to 
determine  accurate  estimates  of  spectral  indices.  We  have  presented  both 
analytical  arguments  and  numerical  experiments,  which  use  least  squares 
procedures  for  simple  geometric  representations  of  scan  functions.  For 
two-region  figures,  such  as  the  top  hat  or  circular  arc,  we  find  that  an 
accurate  estimate  of  the  spectral  index  can  be  obtained  with  256  or  larger 
mode  numbers,  over  a  dynamical  range  of  30-50  db.  However,  for  three- 
spectral  region  figures,  even  with  a  moderate  separation  of  scales,  an 
excess  of  2048  modes  is  required  to  obtain  good  estimates  of  the  spectral 
indices.  Here,  a  typical  dynamical  range  is  80  db.  We  have  Investigated 
various  sources  of  errors  including,  for  example,  the  number  of  modes  in 
the  data  set,  aliasing,  and  transitional  nulls.  It  is  possible  that 
fitting  functions,  which  are  based  on  physical  considerations,  will  reduce 
the  data  set  required  to  obtain  accurate  estimates  of  spectral  indices. 

There  is  a  practical  lesson  to  be  learned  from  this  study.  Spectral 
Indices  determined  from  data  cannot  safely  be  used  beyond  the  range  in  k. 
for  which  they  have  been  measured.  We  have  shown  that  even  under  well 
controlled  conditions  large  errors  in  spectral  indices  can  occur.  To 
extrapolate  these  Indices  beyond  the  dynamical  range  for  which  they  were 
obtained  can  produce  power  level  predictions  that  are  incorrect  by  orders 
of  magnitude.  From  a  systems  perspective  one  could  conceivably  take  a  very 
conservative  point  of  view,  a  worst  case  hypothesis,  and  assume  that 
wherever  the  atmosphere  structures,  the  spatial  power  spectrum  of 
irregularities  falls  off  like  k  .  However,  this  would  be  unduly 
constraining  to  Instrument  designers.  The  fact  is  that  we  do  know  a  fair 
amount  about  the  likely  shapes  of  striatlons  under  HANE  conditions.  We 


second  maximum 


drawn  through  the  relevant  second  maxina 


APPENDIX  A.  TRANSFORM  OF  A  SPECIAL  FUNCTION 


We  consider  Che  function 


f(x)  -  1 


Ixl  <  a, 


f(x)  -  (1  -  -  a)^)^. 


a  <  1x1  <  b, 


f(x)  -  0 


Ixl  >  b, 


where  A  ~  b  '  a  and  t  >  1.  Since  it  is  symmetric  the  transform  can  be 
written  as 


f(k)  -  I  +  I*  +  2k  slo  ka, 


where  the  last  term  is  associated  with  the  region  Ixl  <  a  and 


I  .  .-‘'“I  X 


t+1/2' 


The  last  is  obtained  from  [5]  Sec.  4*3,  Eq.  12  and  the  use  of 
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where  J_  and  I_  are  the  usual  Bessel  functions  and  L  and  H  are  the  Struve 

la  al 

functions. 

It  can  be  shown  that  If  A  »  (b+a)/2,  |f(k)|  has  two  well-separated 
spectral  regions  where  the  spectral  Index  of  the  envelope  Is  pj^  »  2  and  p2 
■  2(t+1),  respectively.  Thus,  If  t  »  -  1  +  e  then  P2  “  2e,  nearly  a  flat 


spectrum. 
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F.  FA JEN 

OICY 

ATTN 

M.  SCHEI3E 

OICY 

ATTN 

CONRAD  L.  LONGMIRE 

OICY 

ATTN 

B.  WHITE 

MISSION  RESEARCH  CORP. 

1720  RANDOLPH  ROAD,  S.E. 
ALBUQUERQUE,  NEW  MEXICO  37106 
OICY  R.  STELLIMGWERl 
OICY  M.  ALME 
OICY  L.  WRIGHT 

MITRE  CORPORATION,  THE 
F.O.  BOX  208 
BEDFORD,  MA  01730 

OICY  ATTN  JOHN  MORGANSTEPJl 
OICY  ATTN  G.  HARDING 
OICY  ATTN  C.E.  CALLAHAN 

MITRE  CORP 

WESTCATE  RESEARCH  PARK 
1820  DOLLY  MADISON  BLVD 
MCLEAN,  VA  22101 


PENNSYLVANIA  STATE  UNIVERSITY 
IONOSPHERE  RESEARCH  LAB 
318  ELECTRICAL  ENGINEERING  EAST 
UNIVERSITY  PARK,  PA  16802 
(NO  CLASS  TO  THIS  ADDRESS) 

OICY  ATTN  IONOSPHERIC  RESEARCH  LAB 

PHOTOHETRICS .  INC. 

4  ARROW  DRIVE 
WOBURN,  MA  01801 

OICY  ATTN  IRVING  L.  KOFSKY 

PHYSICAL  DYNAMICS,  INC. 

P.O.  BOX  3027 
BELLEVUE,  UA  98009 

OICY  ATTN  E.J.  FREMOUH 

PHYSICAL  DYNAMICS,  INC. 

P.O.  BOX  10367 
OAKLAND,  CA  94610 
'  ATTN  A.  THOMSON 

R  i  D  ASSOCIATES 

P.O.  BOX  9695 

MARINA  DEL  REY,  CA  90291 


OICY 

ATTN 

FORREST  GILMORE 

OICY 

ATTN 

WILLI.VM  3.  WRIGHT, 

OICY 

ATTN 

ROBERT  F.  LELEVTER 

OICY 

ATTN 

WILLIAM  J.  KARZAS 

OICY 

ATTN 

H.  ORY 

OICY 

ATTN 

C.  MACDONALD 

OICY 

ATTN 

R.  TURCO 

OICY 

ATTN 

L.  DeR,\ND 

CICY 

ATTN 

W.  TSAI 

RASH)  CORPORATION,  THE 
1700  MAIN  STREET 
SANTA  MONICA,  CA  90406 
OICY  ATTN  CULLEN  CRAIN 
OICY  ATTN  ED  3EDR0ZIAN 

RAYTHEON  CO. 

528  BOSTON  POST  ROAD 


OICY  ATTN  W.  HALL 
OICY  ATTN  W.  FOSTER 


SUDBURY,  MA  01776 

OICY  ATTN  BARBARA  ADAMS 


PACIFIC-SIERRA  RESEARCH  CORP 
12340  SANTA  MONICA  BLVD. 

LOS  ANGELES,  C^  90025 

OICY  ATTN  E.C.  FIELD,  JR. 


RIVERSIDE  RESEARCH  INSTITUTE 
330  WEST  42nd  STREET 
NEW  YORK,  NY  10036 

OICY  ATTN  VINCE  TRAPANI 


37 


* 


SCIENCE  APPLICATIONS,  ISC. 
1150  PROSPECT  PLAZA 
LA  JOLLA,  CA  92037 

OICY  ATTN  LEWIS  M.  LINSON 
OICY  ATTN  DANIEL  A.  HAMLIN 
OICY  ATTN  E.  FRIEMAN 
OICY  ATTN  E.A.  STRAKER 
OICY  ATTN  CURTIS  A.  SMITH 
OICY  ATTN  JACK  MCDOUGALL 

SCIENCE  APPLICATIONS,  INC 
1710  GOODRI0GE  DR. 

MCLEAN,  VA  22102 
ATTN:  J.  COCKAYNE 

SRI  INTERNATIONAL 
333  RATCNSWOOD  AVENUE 
MENXO  PARK,  CA  9A025 


OICY 

ATTN 

DONALD  llEILSON 

OICY 

ATTN 

ALAN  BURNS 

OICY 

ATTN 

G.  SMITH 

OICY 

ATTN 

R.  TSUMODA 

OICY 

ATTN 

DAVID  A.  JOHNSON 

OICY 

ATTN 

WALTER  G.  CHESiniT 

OICY 

ATTN 

CHARLES  L.  RINO 

OICY 

ATTN 

WALTER  JAYE 

OICY 

ATTN 

J.  VICKREY 

OICY 

ATTN 

RAY  L.  LEADABR.AND 

OICY 

ATTN 

G.  CARPENTER 

OICY 

ATTN 

G.  PRICE 

OICY 

ATTN 

R.  LIiaNGSTON 

OICY 

ATTN 

V.  GONZALES 

OICY 

ATTN 

D.  MC0A.N1EL 

TECHNOLOGY  INTERNATIONAL  CORP 
75  WIGGINS  AVENUE 
BEDFORD,  >IA  01730 

OICY  ATTN  H.P.  BOQUIST 

TOYON  RESEARCH  CO. 

P.O.  Box  689(T 
SA-'flA  BARBARA,  CA  93111 
OICY  ATTN  JOHN  ISE,  JR. 
OICY  ATTN  JOEL  GARB.ARINO 


TRW  DEFENSE  &  SPACE  SYS  GROUP 
ONE  SPACE  PARK 
REDONDO  BEACH,  CA  90278 
OICY  .  ATTN  R.  K.  PLEBUCH 
OICY  ATTN  S.  ALTSCHULER 
OICY  ATTN  D.  DEE 
OICY  ATTN  D/  STOCKWELL 
SNTF/1575 

VISIDYNE 

SOUTH  BEDFORD  STREET 
BURLINGTON,  MASS  01803 
OICY  ATTN  W.  REIDY 
OICY  ATTN  J.  CARPENTER 
OICY  ATTN  C.  HUMPHREY 


